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Assemblies of purely repulsive and frictionless particles, display very curious elastic properties 
near the jamming transition where the overlap between particles vanish. Although such systems 
do not contain the long and cross-linked polymeric chains characterizing a rubber, they behave 
macroscopically in a similar way: the shear modulus G can become negligible in comparison with 
the bulk modulus B, as observed in emulsions^. Numerics^ have established a scaling relationship 
between the ratio of elastic moduli and the compression in systems of athermal compressed elastic 
particles. After reviewing recent theoretical results on the microscopic structure of such packings 
(in particular the coordination number z, the average number of contact per particle) I will propose 
an explanation for these observations, and explain why the arguments apply both to soft athermal 
i— H , particles and to colloids where elasticity has an entropic nature. 

hj : 

"sj - \ I. INTRODUCTION 

i— i. Crystalline lattices are invariant under translation, implying that vibrational modes are plane waves. From those 
excitations, one can build a theory of energy transport and a theory of elasticity. In amorphous solids, this symmetry 
breaks down. Although at large length scales a continuous (and translationary invariant) description is a good 
approximation, this fails at small length scales, where the disorder has strong effects. At these scales, various properties 
of amorphous solids, such as energy transport, low- frequency excitations in glasses or force propagation in granular 

S matter, are not yet satisfyingly understood, and are active fields of research. One inherent difficulty in the study of 
these phenomena is that the length scales at play are typically moderate, of the order of ten particle sizes or less. 
'■^J . This makes it harder to test and distinguish clearly the consequences of different theories. Finding a system where 
£h ' length scales can be large and controlled may therefore be extremely useful. 

Emulsion experiments^, followed with theoretical arguments2ii£>£ suggested that the "jamming transition" where 
i ^ i ' repulsive, short-range particles are just in contact corresponds to a critical point. This idea was latter substantiated by 
the findings that the elastic moduli 2 , the vibrational spectrum 2 , the microscopic structure 2 ^ and force propagation- 
C^l ' display scaling behavior near the jamming threshold. At that point, although the system is amorphous and isotropic, 
it cannot be described as a continuous elastic body on any length scaled. Because the strong effects of disorder 
occurs already at large length scales near this critical point, this model system is a lens allowing to probe in a 
stringent manner the properties of amorphous solids, and the effects of disorder. It has enabled to build and tests^^ 
^-j- , a theory for the excess low-vibrational modes (the so-called Boson Peak) found in amorphous solids, which applies 
as well to covalent glasses^ and to model systems of attractive glasse o 10 ' 11 . This line of thought also permits to 
relate microscopic structure and some aspects of the dynamics near the glass transition of hard spheres^ 2 -, despite 
the fact that the glass transition occurs empirically at a packing fraction significantly smaller than those of jammed 
! configurations. 

In what follows I will focus on the elastic moduli, as micro-gels are potentially a good system to vary the packing 
fraction around the jamming threshold to test predictions on elasticity. Near the jamming threshold it is found 
numerically that both for particles interacting with a harmonic potential (used to model emulsions) or a Hertzian 
potential (used to model elastic particles)^, the system is "almost" a liquid, its shear modulus becomes negligible in 
^ , comparison with the bulk modulus: 

% ~ ip/B) 1 ' 2 (1) 

where p is the pressure. In a gel, two different phenomena govern the shear and the bulk modulus: the elasticity of 
the polymeric network and the compressibility of the solvent, respectively. This can cause the two elastic moduli to 
be very different. Nevertheless this behavior is unusual for a solid made of identical particles interacting with a radial 
interaction, where shear and bulk modulus are both induced by local interactions, and are generally comparable in 
amplitude. The theoretical argument, presented in Section llll| will explain this observation and unravel a peculiar 
aspect of the elasticity near the jamming threshold: the elastic moduli can depend enormously on the stress applied 
on the system before the response is measured, the so-called pre-stress. Before discussing this argument, I will start 
by reviewing recent results on the geometry of the packings near the jamming threshold. 
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II. STRUCTURE AND MECHANICAL STABILITY 

Amorphous solids are typically out-of-equilibrium, and their structures a priori depend on their history. It may thus 
seem hard to infer their structure without a detailed description of the way they were made. There is, nevertheless, a 
limiting case that turns out to be conceptually important: when these systems are prepared via a very rapid quench 
from a fluid phase. In this situation, we expect that as soon as the liquid finds some meta-stable states, it remains in 
those states. For slower quenches the system will depart from such states, but there exists evidences that this effect 
is small, at least for hard particles and for quenches as slow as what is typically achieved numerically 12 -. Thus we are 
particularly interested in configurations that are marginally stable, i.e mechanically stable, but which are very close 
to yield. What are those configurations for an assembly of elastic particles? 

Maxwell^ studied mechanical stability in the context of engineering structures, and he found out that the key 
parameter is the coordination z, the average number of interactions per particle. For example, for a network of point 
particle connected via springs, he showed that stability requires z > z c = 2d, where d is the spatial dimension of the 
system. His argument goes as follows: consider a set of N points interacting with N c springs at rest of stiffness k. 
The expansion for the energy may be written: 

SE = J2 |[(^ " ' *kj] 2 + °( SR2 ) (2) 

(ij) 

where the sum is made over all springs, is the unit vector going from i to j, and SRi is the displacement of particle 
i. A system is floppy, i.e. not mechanically stable, if it can be deformed without energy cost, that is if there is a 
displacement field for which SE — 0, or equivalently (SRi — SRj) ■ ~ Vij. If the spatial dimension is d, this linear 
system has Nd degrees of freedom and N c = Nz/2 equations, and therefore there are always non-trivial solutions if 
Nd > N c , that is if z < 2d = z c . To be mechanically stable a network must therefore have z > 2d. 

It turns out that under compression the criterion of rigidity becomes more demanding. Here we shall derive this 
criterion in a simple model, as it yields the correct and more general result. The derivation for a generic amorphous 
packing is made ir>i^. Consider a square lattice made of springs of rest length ^o, which defines our unit length. It just 
satisfies the Maxwell criterion, since z = 4. Now, we add randomly a density Sz of springs connecting second neighbors, 
represented in blue in Fig([T]), such that the coordination is z = z c + Sz. We add them in a rather homogeneous manner, 
so that there are not large regions without blue springs. The typical distance between two blue springs in a given row 
or column is of order Iq/Sz. Dividing this length by the mesh size lg we define the dimensionless number I* = l/5z . 

How much pressure can this system sustain before collapsing? For a system to be mechanically stable, all collective 
displacements need to have a positive energetic cost. It turns out that the first modes to collapse are of the type of 
the red displacement mode represented in Fig((T]): they correspond to the longitudinal mode of wavelength IqI* of a 
segment of springs contained between two blue, diagonal springs. These modes have a displacement field of the form 
SRi = 2Xsva(Tri/l*)/\/T*e x , where i labels the particles along a segment and runs between and I*, e x is the unit 
vector in the direction of the line, and X is the amplitude of the mode, X = 1 for a normalized mode. In the absence 
of pressure p, the energy of this mode comes only from the springs of the segment. In the limit of large /* a Taylor 
expansion of the displacement in Eq. ([2]) gives the well-know result for the energy of the lowest-frequency mode of a 
line of springs SE ~ kX 2 /l* 2 . When p > 0, all the springs now carry a force / ~ pig. The energy expansion contains 
other terms not indicated in Eq.((2]^, whose effect can be estimated quantitatively as follows. When particles are 
displaced along a longitudinal mode such as the one represented in Fig.|T]), the force of each spring directly connected 
and transverse to the segment considered ( see Fig([T]) ) now produces a work equal to / times the elongation of the 
spring. This elongation is simply SR 2 /Iq following Pythagoras' theorem. Summing on all the springs transverse to the 
segment leads to a work of order fX 2 /lo w pX 2 . This gives finally for the energy of the mode SE ~ kX 2 /I* 2 — pX 2 , 
where numerical pre-factors are omitted. Stability requires SE > 0, implying that k/l* 2 > p, or Sz > (p/fc) 1 / 2 . As we 
shall see below, for repulsive particles the bulk modulus B simply follows B ~ fc, and we can rewrite our result as: 

Sz > (p/B) 1 ' 2 (3) 

Physically, this results signifies that pressure has a destabilizing effect, which needs to be counterbalanced by the 
creation of more contacts to maintain elastic stability. The result is more general and for sphere packing one also 
gets Sz = z — z c > (p/B) 1 / 2 . p/B is a measure of the contact strain, and for all interaction potential of interest 
(approximatively harmonic for emulsion or hertzian for elastic body), one finds that p/B ~ (<f> — (j)c), where <f) c is the 
packing fraction at the jamming threshold (this comes from the fact that B = dp/d(j>, leading to p/B ~ (<p — </> c ) if 
power-law behavior are assumed), so that Eq.© can be rewritten Sz > (0— 0c) 1 / 2 . Packings generated numerically^ 
are consistent with the equality of this inequality, supporting that the packing lie indeed close to marginal stability: 
these systems have just enough contacts to counterbalance the destabilizing effect of compression. 
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FIG. 1: Square lattice with a density per particle Sz of additional diagonal springs, represented in blue. I* ~ 1/Sz is the typical 
(dimensionless) distance of the segments contained between two diagonal springs on a given row or column. The red arrows 
represent the longitudinal mode of wavelength ~ I* of such a segment: SRi oc sm(m/l*)e x , following the notation introduced in 
the text. The dashed line exemplifies the deformation of a spring transverse and directly connected to the segment considered, 
it is elongated by the longitudinal vibration of this segment. When the pressure is positive and contacts are under compression, 
this elongation lowers the energy those springs contain. This leads to an instability when 5z becomes smaller than a quantity 
proportional to the square root of the contact strain, of order p/B. 

III. ELASTIC MODULI 

Thus, an assembly of repulsive particles close to the jamming threshold presents a vanishing excess-coordination 
with respect to the jamming threshold: Sz — > as the strain p/B — > 0. Why the proximity of the Maxwell bound 
should impose that shearing becomes much softer than compressing is not obvious a priori. Our argument will show 
that this behavior is specific to purely repulsive systems, and that it is in fact possible to create spring networks where 
the behavior is opposite, for which B/G — > 0. More generally, we shall see that weakly-connected systems (i.e. z close 
to z c ) are soft to most imposed strains, for which they display a dimensionless stiffness or elastic constant of order Sz 
(or zero if the system is floppy, i.e. Sz < 0), with one exception: if the infinitesimal strain is imposed in the direction 
of the stress sustained by the system in its reference configuration, the "pre-stress" . In this case the response of the 
system is not soft and the elastic constant found is similar to the one of a well-connected solid. A simple example is 
given by the zig-zag chain of springs of Fig([2j . When no forces are applied, the chain is in a zig-zag configuration. 
When pulled by the two ends, it yields freely up to the point where the line becomes straight and contact forces 
appear in the springs. At that point, the system does not yield freely anymore, it is stiff. In what follows we shall 
show how this simple idea applies to more complex networks. 




FIG. 2: Zig-zag chain of springs. When pulled at the tips, such a chain yields without stiffness, until the line is straight 
and forces appears in the contacts. For larger strain the stiffness jumps to a finite value. More generally weakly-coordinated 
networks (with a moderate Sz) are soft when an infinitesimal strain is imposed in a generic direction, but are stiff if this 
direction corresponds to the stress (called sometimes pre-stress) sustained by the system. 



A. Force balance and contact deformation operators 

We are left with determining the elastic moduli of packing of a given coordination. We need to introduce some 
amount of linear algebra. The argument proposed is a shortened version of Ref. 10 . Many of the properties discussed 
below concern the response of the system to external forces. It proves convenient to consider our system under the 
influence of an arbitrary set of forces fj acting on all particles i. At equilibrium forces are balanced on each particle 
i: 

^2fijfta=Fi (4) 

U) 
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where /y is the compression in the contact < ij >, the sum is on all the particles j in contact with particle i, and Fi 
the external force applied on i. is the unit vector going from i to j. This linear equation can be written: 

T|f) = |F) (5) 

|f) is the vector of contact tensions and has N c components. F is the vector of external forces. Its dimension is^ Nd, 
since there are d degrees of freedom for the external force on each particle. Therefore T is an Nd x N c matrix. In 
what follows our notations shall be lower-case for the contact space of dimension N c , and upper-case for the particles 
positions space of dimension Nd. 

Another important linear operator describes the change of distances between particles for a given set of displacements 

(SRi - SRj).^ = Srij (6) 

It can be written as: 



S\SR) = \6r) (7) 

where Sr = {Sr^} is the set of distance changes for all contacts. S is an Nd by N c matrix^. 

As was shown e.g. T and S are transpose of each other, as we now prove. At equilibrium, any force field 

applied should cost no energy at first order. This is the virtual work theorem, which reads here: 

J2 SRi ■ Fi - J2 Sr ijfij = m\F) - (<5r|f) = (8) 

i ij 

where we used the scalar product notation notation (<5R|F) = Yli^Ri " F% an d (^ r |f) = Ylij S r ijfij- Applying the 
definitions of S and T in © we get: 

(t\S\SR) = (SR\T\t) (9) 
Introducing the transpose notation this is equivalent to: 

s = r t (io) 



B. Energy expansion and virtual force field 



For simplicity of notation we shall think about rather mono-disperse emulsions, where the interaction potential is 
approximatively harmonic, of stiffness k. The energy expansion follows Eq.(2)2^, which can be rewritten: 

5E = kJ2^ j = ^{SSK\SSR) (11) 

<ij> 

Our system is equivalent to a set of point particles interacting with springs. In order to study its elasticity, it turns 
out to be convenient to consider the responses that follow arbitrary changes of rest length of these springs. This is 
in fact equivalent to imposing dipoles of force. As we shall see, the response to shear or compression can also be 
easily expressed in terms of changes of rest length. We impose an infinitesimal change of rest length on every contact 
y = {uij}. The energy and the displacement field are given by the minimization of: 

k k 
SE= - mmidru -ya) 2 = - mm(S6R.-y\S8~R-y) (12) 

Obviously if S was spanning its image space, we would have SE — : one could always find a displacement |<5R) that 
leads to a change of distances between particles in contact exactly equal to |y). As we said, S is a N c x Nd matrix. If 
N c < Nd, S indeed spans its image space, and the energy associated with any strain |y) is zero: the system is floppy. 
In the other case, if N c > Nd, there are N c — Nd = NSz/2 relations of dependency among the columns of S. One 
can choose a basis of NSz/2 vectors |a p ), with 1 < p < Ndz/2, in the space of \6r) such that: 



(a p |5 = 



(13) 
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The |a p ) are orthogonal to all the vectors <S|<5R), for any displacement field \SR). Transposing this relation we have: 

T\aP) = (14) 

which indicates that all the vectors in the space of the |a p ) satisfy force balance without external force but no 
others. The |a p ) live in the contact-force space, and henceforth we shall denote them |a p ) = |f p ) = {ffj}- In the 
following we consider an orthogonal unit basis: 

(F|fp') = ^/P^'= V (15) 

y 

We can decompose any |y) as: 

y = y^+ E ( 16 ) 

p=l...N% 

y , the part of |y) orthogonal to the |f p ), is spanned by the matrix S, and therefore does not contribute to the energy 
when the minimization of Eq. (|12[) is performed. In other words, there exists a displacement field <5Ro which leads to 
a strain y 1 - = SSRq. On the other hand, the strain field corresponding to the second term in Eq. fTI)]) is orthogonal 
to the space generated by S, and cannot be accommodated by any displacements. The minimum of Eq. (j!2[) therefore 
occurs in <5Ro and one finds for the energy: 

5E= k - ]T (F|y) 2 (17) 

We now discuss properties of the contact force fields that we shall use to derive the scaling of the bulk and the 
shear moduli. Only one vector of the vector space of the force fields |f p ) solutions of Eq.© without external force 
is the real set of contact forces that supports the system, and that could be observed empirically. This vector is 
denoted If 1 ). The rest of the basis |f p ) with 1 are also solutions of Eq. (O without external force. Nevertheless, 
there are not "physical" solutions for the interaction potential chosen. Thus we shall call them virtual. If 1 ) verifies 
the following properties: (i) In a system with repulsive interaction, as we consider here, all the contact forces are 
compressive and therefore > for all contacts, (ii) It is well known from simulations and experiments that the 
distribution of contact forces is roughly exponential, or compressed exponential (see for example 2 for simulations in 
the frictionless case). This implies that the fluctuations of the contact forces are of order of the average value, leading 
t° (flj) 2 ~ ((fij) 2 ) — for a normalized force field. Thus we may introduce a constant cq such that: 

<4>= c o-±= (18) 

Now we turn to the properties of the virtual forces |f p ): (i) There are no physical constraints on the sign of the 
contacts force for the virtual vectors. Furthermore, the |f p ) must be orthogonal to If 1 ), whose signs of contact forces 
are strictly positive, and where the fluctuations in the contact compression is small. Therefore the virtual force fields 
have roughly as many compressive as tensile contacts. 



C. Elastic moduli 



In our framework it turns out to be convenient to study the response to shear or compression as there are generally 
implemented in simulations. When periodic boundaries conditions are used, an affine strain is first imposed on the 
system. Then the particles are let to relax. In general the affine strain is obtained by changing the boundary condition. 
Consider a 2-dimensional system with periodic boundary: it is a torus. For example a shear strain can be implemented 
by increasing one of the principal radii of the torus and decreasing the other. Then the distance between particles in 
contact increases or decreases depending on the direction of the contact. In fact, this procedure of change of boundary 
conditions is formally equivalent to a change of the metric of the system. If the metric is changed from identity / to 
the constant metric G = I + U + lP , where U is the imposed infinitesimal global strain, the length of a vector 51q 
becomes SI, such that SI 2 — 6l ■ G ■ 6l - Using this expression with SIq — Rij = Rj — Ri one deduces the change of 
distance between two particles: 



(19) 
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Near jamming, for mono-disperse particle of diameter Iq, Rij ~ loftij and therefore Sr^ ~ Zo^ij ' {U + U ) • Hij/2. 
Formally, such a change of metric is strictly equivalent to a change of the rest length of the springs with j/y = br^ . 
Incidentally Eq. (fl7| can be used to compute the energy of such strain. 

For a compression (U + U t )/2 = —el where I is the identity matrix and e is the magnitude of the strain. Eq. (fT7|) 
becomes: 

6 -§ = \(E-^? + \ E (-Ei£) 2 (20) 

In the first sum all the terms have the same sign for a purely repulsive system, and this term leads to the strongest 
contribution. We have: 

§ > (£ -e#) a = e 2 {N c ( flJ )) 2 = e 2 c 2 N c (21) 

ij 

On the other hand, 8E is certainly smaller than an affine compression whose energy also goes as ke 2 N. Therefore we 
find that: 

5E ~ kNe 2 l 2 (22) 
B = ~ kit 2 (23) 

which does not depend on Sz. Here V ~ NIq is the volume of the system. The bulk modulus of an harmonic system 
jumps from in the fluid phase toward a constant when the system becomes jammed, as observed in the simulations. 
From this result follows that p ~ B(<fi—(j) c ) ~ (cf) — (f> c ). Note that this result holds only for purely repulsive systems. If 
the potential has an attractive component, and if the pressure is set to zero, the real force field If 1 ) presents as many 
negative contact forces as positive ones, and this term does not lead to a particularly large contribution in Eq. (|2"0|) . 
In this case, one expects generically to recover for the bulk modulus the result valid for the shear modulus, that we 
derive in the next section. 

If a pure a shear strain is imposed, the tensor (U + U t )/2 is traceless. Let be e the largest eigenvalue (in absolute 
value). The change of distance of two particles in contact due to a pure shear Srij is a number of zero average if the 
system is isotropic, and fluctuates between +e and — e depending on the orientation of Rij. Eq. (|17[) becomes: 

* E =\ E (E/£^) 2 ( 24 ) 

p=l,..«|£ Ij 

Each term in the summation gives on average: 

((E f?i Sr vf) - E((/5) 2( H> + E (f!if P m Jn 3 Sr mn ) (25) 

ij ij mn^ij 

= E<(/5) 2 ><*4>+ E (/5/£«^Amn> (26) 
ij mn^ij 

For a generic imposed strain, Srij and Sr mn are not correlated, so that the terms (ffj fmn^ r ij^ r mn) vanish in Eq|25l 
For a pure shear, spatial correlations between Srij and 5r mn exist, and in order to estimate the shear modulus one 
must make an extra assumption on the nature of the disorder of the network being considered. For isotropic random 
sphere packings one expects the contact network to be strongly disordered and the propagation of forces to be strongly 
scattered. We shall therefore assume that the contact forces fields only present weak spatial correlations, as is indeed 
observed numerically^ (note that it is possible to build networks with reduced disorder where this assumption breaks 
down, for example in the example shown in Fig.{l} where long straight lines are present in the microscopic structure. 
If this system is elongated in the direction of these lines, the restoring force will be large independently of the excess 
coordination Sz). Following this assumption for packing of particles, we shall neglect the terms (ffj fmn^ r ij^ r mn) 
when mn ^ ij. Concerning the diagonal terms, one has Srfj « l^e 2 while J2(ffj) 2 = 1 by construction. Thus each 
term in the p summation is of order IqC 2 ■ 1, and: 

5E - kN8ze 2 l 2 (27) 
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Note that this estimation for the energy will apply for a generic imposed contact strain field |y), since in general such 
a field would only project weakly on the real force field If 1 ). Finally Eq (j2"7"|) implies: 

SE 

G = ~ kSzl 2 - d (28) 

Eqs. (|28l22p lead to the result G/B ~ Sz. In the context of frictionless particles, one may use Eq.([3]) to recover the 
observed result G/B ~ (p/B) 1 / 2 . 



IV. SUMMARY AND CONCLUSION 



We have shown that weakly-coordinated elastic networks display an exotic elastic behavior: they are very soft 
and present an elastic constant proportional to Sz (of order 8zkl^~ 2 , rather than kl^~ 2 expected for a usual solid) 
for generic deformations, but recover a usual, normally-connected behavior when the strain is imposed along the 
pre-stress. These results apply to interacting particles, with the difference that the typical stiffness k now depends on 
compression. For emulsion the potential is nearly harmonic and the dependence of k with compression is negligible. 
For elastic particles interacting with a Herztian potential, one gets k ~ p 1 ^ 3 . For a hard sphere glass, which lies below 
jamming <j) < <j> c , the interaction is purely entropic, and is studied in ReiM. One finds k/ksT ~ (p/fc^T) 2 . 

A consequence of our analysis is that systems made of repulsive particles, which are always under pressure, are stiff 
to compression, with respect to the shear modulus that can be tiny if the contact strain is small. At a qualitative level 
this has been observed empirically in emulsions^, glass beadsi£ and sancUI. Nevertheless as far as scaling exponents 
are concerned this results has been backed only numerically, without 1 -! 2 " and with frictio n 16 i 18 . In order to observe this 
phenomenon quantitatively in nature, emulsions and micro-gel are presumably good candidates. Ideally one would 
like to have a system where (i) friction is small, as it limits how close one can get from criticality-i^ (ii) thermal effects 
are small near the jamming threshold, for the same reason. This requires to consider sufficiently large particles, (iii) 
The osmotic pressure is controlled to vary the distance to threshold reliably. 

Such an experimental system may enable to address additional intriguing questions, in particular how amorphous 
solids yield under shear. From our analysis we predict that as a shear strain is imposed, and as the shear stress 
builds up, the structure must stiffen. For example if the shear stress can be made of the order of the compression 
near the jamming threshold -which remains to be established- the shear modulus must become of the order of the 
bulk modulus. Equivalently, as force chains orientate to hold the shear stress, the system becomes much stiffer in the 
directions where force chains align. At a qualitative level such a stiffening has been observed numerically 2 ^ with elastic 
particles, and we have argued that this effect occurs more generally in various systems such as in gels of semi-flexible 
polymers 2 ^. The small value of the shear modulus reflects the presence of soft collective modes in the amorphous 
structure 2 ^, which stop coupling to the applied strain when stiffening occurs. It would be interesting to investigate if 
these modes, and the stiffening they generate, are causally related to the yielding of the structure. The scaling of the 
maximal strain with compression before yielding would be very informative as well to address this issue. 
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More precisely, its dimension is (Nd — d(d + l)/2). The term d(d + l)/2 corresponds to the constraints on the total torques 
and forces that must be zero at equilibrium. 

If one removes the global translations or rotations from the displacement fields, which obviously do not change any inter- 
particle distance, 5 becomes a (Nd — d(d + l)/2) by N c matrix. 

When forces are present in the contact, as is the case at finite pressure, another term enters the energy expansion, see^. It 
does not affect the scaling of the elastic moduli, and we shall ignore it. 



